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ABSTRACT 

We present light curves and periodograms for 27 stars in the young Upper Scorpius 
association (age=ll± IMyr) obtained with the Kepler spacecraft. This association is 
only the second stellar grouping to host several pulsating pre-main sequence (PMS) 
stars which have been observed from space. 

From an analysis of the periodograms, we identify six 6 Scuti variables and one 
7 Doradus star. These are most likely PMS stars or else very close to the zero-age 
main sequence. Four of the <5 Scuti variables were observed in short-cadence mode, 
which allows us to resolve the entire frequency spectrum. For these four stars, we are 
able to infer some qualitative information concerning their ages. For the remaining 
two 6 Scuti stars, only long-cadence data are available, which means that some of the 
frequencies are likely to be aliases. One of the stars appears to be a rotational variable 
in a hierarchical triple system. This is a particularly important object, as it allows 
the possibility of an accurate mass determination when radial velocity observations 
become available. We also report on new high-resolution echelle spectra obtained for 
some of the stars of our sample. 

Key words: Stars pre-main sequence - Stars: oscillations - Stars: variables: delta 
Scuti - Stars: variables: Herbig Ae/Be - Stars: fundamental parameters 


1 INTRODUCTION 

A complete understanding of the star formation process re¬ 
quires the ability to predict how the properties of young stars 
depend on their initial conditions. Despite significant efforts, 
a key observational problem remains: the reliable determi¬ 
nation of the masses and ages of the youngest pre-main- 
sequence (PMS) stars once they become optically visible, 
after the main accretion phase. Stellar masses are needed 
to investigate the shape of the initial mass function (IMF) 
and its possible dependence on the properties of the parent 
cloud. However, obtaining reliable stellar masses is a very 
difficult task. The only direct method is to use eclipsing bi¬ 
nary systems where individual masses of the components 
can be obtained with great accuracy from the radial ve¬ 
locity and light curve. However, eclipsing systems for PMS 
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stars are rare. Another commonly adopted procedure is to 
locate the stars in the Hertzpung-Russell (HR) diagram, us¬ 
ing estimates of their effective temperatures and luminosi¬ 
ties. Masses can be obtained by comparing their position in 
the HR diagram with theoretical evolutionary tracks, assum¬ 
ing the chemical composition (Hillenbrand fc White|[M04 \. 


Recently, considerable advances have been made in our 


understanding of PMS evolution (Palla & Stabler 1999 


Siess, Dufour, fc Forestini|2000||Baraffe et al.|2002||Tognelli,| 


Prada Moroni, fc Degl’InnocentipOll I. However, differences 
in the evofutionary tracks remain due to the different treat¬ 
ment of convection, different opacities and differences in the 
zero-point of the calculated ages. This results in consider¬ 
able differences in mass estimates of single stars from their 
effective temperatures and luminosities. Without adequate 
observational constraints, it is impossible to determine the 
correct assumptions in the evolutionary model calculations. 
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For this reason, it is very important to obtain independent 
determinations of PMS masses and ages. 

In this context, intermediate-mass PMS stars (with 
mass in the range 1-8 Mq, also called Herbig Ae/Be stars) 
are particularly useful. After the seminal work by |Marconi| 
& Palla (19981, who established the locus of the theoreti¬ 


cal instability strip in the HR diagram, it has become clear 


that these stars pulsate as S Set variables (see e.g. Ripepi et 
arconi fc Palla|2004| [Ripepi et al.|2006a|b[ 


al.|2002[[2003l|Ma 


Zwintz |2008| and references therein). More recently, high- 
precision photometry from space missions such as MOST 
and CoRoT showed that intermediate-mass PMS stars can 
also pulsate as hybrid S Sct /7 Dor stars | Ripepi et aLpOlT I 
or pure 7 Dor stars (Zwintz et al. 20131. This is interest¬ 


ing because the occurrence of low frequencies in 6 Set stars 
(i.e. hybrid pulsation) was originally thought to be excep¬ 
tional. However, Kepler observations have shown that all 


S Set stars have low frequencies (Balona 20141. The rea¬ 


son why the hybrid concept arose is due to the fact that 
the low frequencies attain sufficient amplitudes to be visible 
from the ground only in a narrow temperature range. The 
mechanism which drives low frequencies in hot S Set stars 
is currently unknown. 

The additional information provided by pulsation can 
be used to construct detailed models of the varying internal 
structure of these stars as they contract toward the main 
sequence, in order to compare the predicted frequencies to 
the observed ones (Ruoppo et al.|[2007 Di Criscienzo et al. 


2008 Casey et al.|2013 l. The physics of PMS stellar evolu¬ 

tionary models can therefore be tested. Several efforts have 
recentfy been made to anafyze the properties of known PMS 
pulsators (see e.g. Zwintz et af.|20l'4 l. However, the anafysis 
is hampered by severaf factors such as the uneven quality of 
the observational data (ground and space time series data 
have very different properties), the lack of precise spectro¬ 
scopically derived stellar parameters, the possibly different 
intrinsic properties of cluster and field PMS pulsators, and 
difficulties in the treatment of rapidly rotating stars. 

In this context, PMS pulsators belonging to clusters 
or associations are of particular interest, since the member 
stars share the same chemical composition, age and inter¬ 
stellar absorption, greatly simplifying the comparison be¬ 
tween theory and observations. This is the case, for example, 
of NGC2264, a well-studied cluster where 9 PMS pulsators 
have been observed from space ( [Zwintz et al.|20l4 l. The Up¬ 
per Scorpius association (USco) is perhaps of even greater 
interest due to the large number of known A-F members 
(de Zeeuw et al. |1999| [Pecaut, Mamajek, fc Bubar||20l^ 
and the availability of accurate individual distances based 
on Hipparcos parallaxes. Moreover, stellar parameters have 
been estimated for most of the member stars. 

USco is part of the larger Scorpius-Centaurus (Sco- 
Cen) association, a young (~5-10 Myr) and relatively close 


(~120-150pc) region of recent star formation (de Brui- 
1999| de Zeeuw et al. 19991 that includes the Up- 


jne 


per Centaurus-Lupus (UCL) and Lower Centaurus-Crux 
(LCC). USco is the youngest of them and contains many 
intermediate-mass members, making it an ideal target to 
investigate the whole extension of the PMS instability strip. 
Indeed, its member stars are sufhciently bright (~ 7 < U < 
10 mag) to acquire extremely precise photometry both from 
the ground and from space. 


Due to the loss of a second spacecraft reaction wheel, 
the Kepler satellite ended data collection in the original 
Cygnus field after four years of continuous monitoring. By 
pointing near the ecliptic plane, the Kepler spacecraft is 
able to minimize pointing drift, acquiring data with much 
reduced photometric precision. This mode of operation is 
called the K2 projeetj^The USco region falls partially within 
the field of view of the Kepler satellite during the K2-C02 
campaign, becoming a primary target for our purposes. In 
this work we present the results of a search for pulsating 
PMS stars among the USco member stars observed by Ke¬ 
pler. 

The structure of the paper is the following: in Section 
2 we describe the Kepler data and the techniques of data 
reduction and analysis; in Section 3 we present the classifi¬ 
cation of the variables; Section 4 reports a the results of the 
spectroscopic follow up for a sample of interesting templates; 
Section 5 includes a discussion about the evolutionary sta¬ 
tus of the target stars; Section 6 reports a discussion of the 
properties of the pulsating stars. Finally, a brief summary 
closes the paper. 


2 KEPLER OBSERVATIONS 


Kepler observations usually consist of 30-min exposures 
(long cadence or LC mode), but a few stars can be observed 
with 1-min exposures (short cadence or SC mode). Obser¬ 
vations of 16 PMS stars in USco, 7 of which were also ob¬ 
served in SC mode, were obtained by Kepler. The K2-C02 
campaign started on 2014 Aug 23 and ended on 2014 Nov 
10, for a total duration of about 78.7 d. The selected stars 


come from the lists of de Bruijne (19991 and Pecaut, Mama¬ 


jek, & Bubar (20121. The spectral types of the selected stars 


range from A2 to F7 (to embrace the entire S Set and 7 Dor 
instability strips, taking also into account possible errors on 
the spectral type values by 1-2 subclasses or more). We also 
searched the Kepler K2 data archive at MAST (Mikulski 
Archive for Space Telescopes) for additional A-F members 
of USco and found an additional 9 objects (2 in SC mode), 
mainly AO stars, giving a total of 27 USco members. The 
targets and their main properties are listed in Table 


2.1 Data reduction and frequency analysis 


Data reduction was carried out by means of the PyKe soft¬ 


ware (Still & Barclay 20121 and also by our own custom 


software. Simple aperture photometry (SAP) was extracted 
from raw pixel files and then corrected for the instrumental 
signatures present in the data. Spacecraft drift is the main 
component leading to increased photometric scatter com¬ 
pared to the original Aep/er field. Using SAP, the brightness 
of a target changes slowly as the star moves from its original 
location on the CCD until the thrusters are applied to bring 
the star back to its original location. The periodogram there¬ 
fore shows strong peaks at about 4.08 d~^ and its harmonics, 
which corresponds to the frequency at which the thrusters 
are applied. It is possible to obtain a greatly improved light 
curve simply by locating these peaks in the periodogram 


^ See http://keplerscience.arc.nasa.gov/K2/ for details. 
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Figure 1. LC light curves for the 27 objects investigated in this paper. Black and red points show data with quality flag equal or greater 
than zero, respectively. For some objects the non zero flag data was omitted because it was too scant. 
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Table 1. Stellar Parameters For A-type Upper Sco Member. The meaning of the different columns is the following: (1) EPIC number; 
(2) HIP identifiers; (3) mode of observation: LC and SC stands for Long and Short cadences, respectively; (4) variable classification; (5) 
if epZer magnitude; (6) Spectral Type; (7) Interstellar absorption; (8) parallax; (9) logarithm of the effective temperature; (10) logarithm of 
the luminosity in solar units. All the data reported in columns (6-10) are from |Pecaut, Mamajek, fc Bub^ | |2012[ |; (11) Disk classification 
according to|Luhman fc Mamajeki (|2012^: d/e=debris/evolved transitional disk; f=fuU disk. 


EPIC 

(1) 

HIP 

(2) 

Mode 

(3) 

Var. Type 

(4) 

KepMag 

mag 

(5) 

ST 

(6) 

Av- 

mag 

(7) 

TV 

mas 

(8) 

log Tetf 

(9) 

logi/Lg 

(10) 

Disk 

(11) 

202842502 

79733 

LC 

_ 

9.054 

AlmA9-F2 

1.25 ± 0.04 

4.51 ± 0.43 

3.964 ± 0.019 

1.48 ± 0.092 

_ 

202876718 

79878 

LC 

- 

7.183 

AOV 

0.00 ± 0.02 

7.35 ± 0.55 

3.980 ± 0.031 

1.407 ± 0.07 

d/e 

203120347 

80586 

LC 

- 

8.211 

F5IV-V 

O.OOiO.lO 

7.02±0.65 

3.814±0.011 

0.93±0.09 

- 

203399155 

80311 

LC 

- 

8.878 

AIV 

0.93 ± 0.05 

5.51 i 0.49 

3.964 ± 0.019 

1.266 ± 0.08 

- 

203660895 

79097 

LC 

- 

9.155 

F4V 

0.48±0.12 

6.99±0.57 

3.822±0.009 

0.90±0.09 

- 

203712541 

78494 

LC 

- 

7.722 

A2mA7-F2 

0.75 ± 0.06 

7.26 ± 0.57 

3.943 ± 0.016 

1.395 ± 0.077 

- 

203774126 

80799 

LC 

- 

7.900 

A3V 

0.25 ± 0.04 

7.97 i 0.62 

3.932 ± 0.012 

1.078 ± 0.068 

d/e 

203931628 

80196 

LC 

DSCT 

9.030 

AlVn 

2.36 ± 0.18 

5.94 ± 0.56 

3.964 ± 0.019 

1.602 ± 0.102 

d/e 

204054556 

79054 

LC 

GDOR 

9.215 

F3V 

0.35 ± 0.05 

6.24 ± 0.53 

3.827 ± 0.005 

0.78 ± 0.08 

d/e 

204076987 

79643 

LC 

- 

9.393 

F3V 

0.56±0.05 

6.63±0.62 

3.827±0.005 

0.67±0.08 

d/e 

204175508 

77960 

LC 

DSCT 

8.410 

A4IV/V 

0.75 ± 0.02 

8.56 ± 0.71 

3.918 ± 0.012 

0.973 ± 0.093 

- 

204222666 

78099 

LC 

- 

7.790 

AOV 

0.59 ± 0.03 

7.42 ± 0.59 

3.980 ± 0.031 

1.367 ± 0.073 

d/e 

204239132 

80238 

LC 

- 

7.848 

A2.5V 

0.74 ± 0.09 

8.00 ± 0.68 

3.937 ± 0.014 

1.389 ± 0.081 

d/e 

204242194 

79250 

LC 

- 

7.829 

A3III/IV 

0.33 ± 0.05 

9.60 i 0.75 

3.932 ± 0.012 

0.963 ± 0.072 

d/e 

204372172 

80088 

LC/SC 

DSCT 

9.000 

A9V 

0.61 ± 0.03 

6.11 ± 0.66 

3.872 ± 0.010 

0.960 ± 0.095 

d/e 

204399980 

79476 

LC/SC 

DSCT/LPV 

8.822 

A8IVe 

1.02 ± 0.40 

7.56 ± 0.61 

3.875 ± 0.010 

1.006 ± 0.175 

f 

204492384 

79644 

LC 

ROT 

10.024 

F6V 

0.53±0.11 

6.30±0.64 

3.S02±0.007 

0.44±0.10 

- 

204494885 

80130 

LC/SC 

DSCT 

8.727 

A9V 

0.66 ± 0.02 

6.41 ± 0.57 

3.872 ± 0.010 

1.089 ± 0.077 

- 

204506777 

78977 

LC/SC 

EA/ROT 

8.752 

F8V 

0.38±0.08 

7.53±0.64 

3.788±0.007 

0.83±0.08 

- 

204638251 

80059 

LC/SC 

DSCT 

8.745 

A7III/IV 

0.56 ± 0.07 

7.74 ± 0.64 

3.892 ± 0.014 

0.907 ± 0.075 

- 

204662993 

79977 

LC/SC 

- 

9.099 

F3V 

0.36±0.05 

7.79±0.66 

3.827±0.005 

0.63±0.08 

d/e 

204810792 

79083 

LC 

- 

8.226 

F3V 

0.76±0.13 

7.16±0.63 

3.S27±0.005 

1.13±0.09 

- 

204926239 

79606 

LC/SC 

- 

9.039 

F8V 

0.81±0.14 

9.19±0.79 

3.788±0.007 

0.65±0.09 

- 

204966512 

80019 

LC 

- 

8.396 

AOV 

1.03 ± 0.07 

7.81 i 0.68 

3.980 ± 0.031 

1.374 ± 0.084 

d/e 

205002311 

82218 

LC 

ROT? 

9.008 

F3V 

0.30±0.02 

8.39±0.72 

3.827±0.005 

0.55±0.08 

d/e 

205089268 

79156 

LC 

- 

8.804 

AOV 

0.60 ± 0.07 

6.43 ± 0.54 

3.980 ± 0.031 

1.366 ± 0.081 

d/e 

205181377 

79124 

LC 

- 

8.630 

AOV 

0.78 ± 0.05 

6.51 ± 0.53 

3.980 ± 0.031 

1.55 ± 0.078 

- 


and fitting and removing a truncated Fourier series. Extrac¬ 
tion of the light curve from the FITS files using SAP, fitting 
and removing the truncated Fourier series was done using 
custom software. 


A more sophisticated approach to correct the light curve 
was proposed by |Vanderburg fc Johnson (20141. In this tech¬ 
nique the non-uniform pixel response function of the Kepler 
detectors is determined by correlating flux measurements 
with the spacecraft’s pointing and removing the dependence. 
This leads to an improvement over raw SAP photometry by 
factors of 2-5, with noise properties qualitatively similar to 
Kepler targets at the same magnitudes. There is evidence 
that the improvement in photometric precision depends on 
each target’s position in the Kepler field of view, with worst 
precision near the edges of the field. Overall, this technique 
restores the median-attainable photometric precision within 
a factor of two of the original Kepler photometric preci¬ 
sion for targets ranging from 10-15 magnitude in the Kepler 
bandpass. This technique is implemented in the kepsff task 
within the PYKE software suite. This is the method used in 
our analysis. Correction of the periodogram of the raw SAP 
photometry as mentioned above is very useful as a check on 
possible artefacts which may arise in the kepsff task. 


The light curves for the 27 targets with LC data are 
shown in Fig. whereas Fig. displays an enlargement of 
the light curves for the 7 stars with SC data. Periodograms 
and frequency extraction were performed using the Pe- 


riod04 software (Lenz & Breger 2005k. The periodograms 


il 


for these stars are shown in Fig. and 

Problems in the low frequency regime may be expected 
given the image drift. It is therefore important to identify 
common frequencies in the periodograms which could be 
of instrumental origin or artefacts of the reduction process. 
Most of the periodograms obtained from applying the [Van-^ 


derburg & Johnson (20141 technique have a peak at fre¬ 


quency u — 0.51 d~^ which can be identified as an instru¬ 
mental frequency. In many stars the peak at half this fre¬ 
quency, ~0.26d~^, is also present. Stars having this feature 
include EPIC 202842502, 203660895, 203931628, 204076987, 
204494885, 204638251, 204810792 and 205089268. Other ar¬ 
tificial peaks are present and these must be recognized in 
order not to be misled in classifying the star as a variable. 
In Fig. and we have shown with dashed lines the fre¬ 
quencies that are likely of instrumental origin. These are the 
harmonics of = 0.51 d~^ and of 4.08d~^. 


3 CLASSIFICATION OF VARIABLES 

As an essential first step, it is necessary to examine the pe¬ 
riodograms and light curves of all stars in order to identify 
pulsating variables and other interesting objects. The only 
pulsating stars in the observed spectral range that can be 
safely identified are S Set (characterized by the presence of 
one or more high frequencies, i.e. frequencies in excess of 
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Figure 1 - continued 


about 5 d“^). 7 Doradus stars, which have multiple frequen¬ 
cies lower than about 5d~^, are difficult to identify with 
certainty. This arises because of the possibility of spurious 
instrumental frequencies in this range as well as irregular 
variability caused by circumstellar material which might be 


present in these young stars. Furthermore, one cannot ig¬ 
nore the possibility that many stars may have spots which 
have different frequencies due to differential rotation or hnite 
lifetimes. Bearing these difficulties in mind, our attempted 
classifications are presented in Table 
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Figure 1 - continued 


Five stars are definitely 5 Set (DSCT) variables: 
EPIC 203931628, 204175508, 204372172, 204494885, and 
204638251. At a first look EPIC 204399980 seems to be an 
irregular long-period variable. However, a closer inspection 
of the SC light curve and its periodogram (Fig. and 


reveals a very likely 5 Set pulsation with frequencies in the 
range 30-40 d“^ superimposed over irregular variations on 
a longer time scale. This kind of photometric behaviour is 
typical of an Herbig Ae star and is caused by the partial ob- 
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Figure 2. The figure shows selected portion of SC light curves for the seven targets observed in this cadence. 


scuration of the star’s photosphere by circumstellar material 
(see Section]^. 

We identify EPIC 204054556 as a 7 Dor (GDOR) vari¬ 
able. EPIC 204492384 could be GDOR, but it looks more like 
a rotational variable (broad multiple peak at ~1.14d“^ and 


harmonic at ~2.28 d“^). The multiple close frequencies at ~ 
1.14 d~^ can be attributed to differential rotation. Similarly, 
EPIC 205002311 might be a rotational variable because the 
harmonic is present, but this classification is uncertain due 
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Figure 3. Fourier transform for the 27 objects with LC data. The dashed red lines show frequencies related to the satellite motion. 


to the low amplitude and the significant noise in this part 
of the spectrum. 

EPIC 204506777 (HD 144548) is a very interesting star. 


It is an eclipsing binary of Algol type (Kiraga 


20121 with 


two small amplitude eclipses. From the periodogram we find 
a main period P = 1.524(1) d. The light curve at this period 
is roughly sinusoidal and could be interpreted as the rotation 
period. There also seems to be a short regular eclipses at P = 


period P = 1.62780 d and usini = 80 ± 5 km s . There are 1.628(1) d (the same period reported by Kiraga 2012 see also 
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Figure 3 - continued 


Fig.[2f. On top of this are deeper eclipses at JD 2064.32 and 
JD 2124.01 and a secondary double eclipse at JD 2089.94 and 
JD 2090.30. It is not easy to interpret this variability, but 
possibly what we are seeing is a third body which eclipses 
a close binary with period 1.628(1) d. The double eclipse 


is of equal depth, indicating that the close binary has equal 
components. This is a very interesting system which deserves 
further study. 
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Figure 3 - continued 


4 SPECTROSCOPIC FOLLOW-UP 

In the previous section we identified several interesting stars 
worth of further analysis, especially the seven pulsating stars 
and the possible triple system (EPIC 20450677). To fully ex¬ 


ploit the potentialities of the K2 photometry, especially for 
the pulsating stars, it is important to derive high-resolution 
spectroscopic estimates of the stellar parameters and of the 
rotational velocities usini. To this aim we searched the lit- 
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Figure 4. Fourier transform for the 7 objects with SC data. The dashed red lines show frequencies related to the satellite motion. 


erature and the available web archives for high-resolution 
spectra of our targets. We found data for two stars. 

• EPIC 204399980 (HIP 79476 = HD 145718) has been 
studied spectroscopically in some detail by [Carmona et ah] 


(20101. The spectrum displays an inverse P Cygni Ha pro¬ 
file, presumably a result of accretion, confirming the Herbig 
Ae nature of this star. They obtained a projected rotational 
velocity, usini = 100 ± lOkms”^ and a centre-of-mass ve¬ 
locity u,. = 0 ± 3kms“^. 
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• A spectrum of EPIC 204506777 (HIP 78977 = 

HD 146897), obtained in 2006, 20 June with the FEROS 
echelle spectrograph on the 2.2-m MPG telescope (resolu¬ 
tion of 48 000, exposure time 90s) was found in the ESO 
archive. 

Nothing was found for the remaining 6 pulsating stars. 
Hence we decided to to observe all of them using the 
new Catania Astrophysical Observatory Spectropolarimeter 
(CAOS) facility. 

The next section describes these new spectroscopic ob¬ 
servations and data analysis for the seven observed stars, 
(see Table 1^ for a list). We note that both EPIC 204399980 
and EPIC 204506777 were re-observed with CAOS in or¬ 
der to search for duplicity in the former star and to com¬ 
pare the spectroscopic results for the latter. We also ob¬ 
served the two stars with SC data which do not pulsate, i.e. 
EPIC 204662993 and EPIC 204926239. Unfortunately, the 
weather conditions did not allow observations of the pul¬ 
sating stars EPIC 203931628, 204054556 and 204175508. 


4.1 CAOS spectra 


CAOS is a fiber fed, high-resolution spectrograph recently 
installed at the Cassegrain focus of the 91-cm telescope of 
the “M. G. Fracastoro” observing station of the Gatania 
Astrophysical Observatory (Mt. Etna, Italy). The spectra 
were obtained in 2015, June and July. Because the targets 
were very low above the horizon, they could only be observed 
for about one hour. We used an exposure time of 3600 s 
for all observations. The typical signal-to-noise ratio of the 
spectra is in the range 40-80 and the resolution R — 45000 
as measured from the ThAr and telluric lines. For details 
on the spectrograph and the data reduction procedure, see 
[Catanzaro et aL] | |2015[ ). 

The spectrum synthesis approach followed in this pa¬ 
per to derive fundamental stellar parameters is based on 


the Unix port of the ATLAS9-SYNTHE suite of codes (Ku 


rucz fc Avrett [1981 Kurucz|1993a|b Sbordone et al.|2004 |. 
This approach has been successfully used in a number of pa¬ 
pers in the recent literature (see e.g. Niemczura et al. pol^ 


Catanzaro fc Ripepi| 


[Hubrig et al. 2014| Doyle et ah] 


2014[ Paunzen et al.||2014f Catanzaro et al. |2015| just to 

quote some recent work). The website by Fiorella Castell0 
contains a description of all the improvements of the Kurucz 
codes and related updates in the input physics performed in 


the latest years. A comparison between SYNTHE (Kurucz 


& Avrett 1981 ) and SYNTH (Valenti & Piskunov 19961 was 


already performed by Catanzaro, Ripepi, & Bruntt (20131, 


who demonstrated that the results provided by the two 
codes are perfectly consistent. We also note that the new 
version of ATLAS, i.e. ATLAS12 ( Kurucz|1997 1, differs from 
the previous ATLAS9 only in the method adopted to com¬ 
pute the opacity. A comparison between ATLAS9 and AT- 
LAS12 synthesis was performed by e.g. [Catanzaro fc Balona| 
(20121, who found that the two codes give consistent re¬ 


sults if the metallicity of the target star is not far from the 
adopted solar one. 

In this paper, since the targets likely have circumstellar 


^ http://www.oact.inaf.it/castelli/ 


material that could partially fill the Ha line, we have derived 
the effective temperatures minimizing the difference between 
observed and synthetic H/? profiles, using as goodness-of-ht 
parameter the defined as 



where N is the total number of points, lobs and Ljt are the in¬ 
tensities of the observed and computed profiles, respectively, 
and Slobs is the photon noise. Synthetic spectra were gener¬ 
ated in three steps: i) we computed LTE atmospheric models 
using the ATLAS9 code; ii) the stellar spectra were synthe¬ 
sized using SYNTHE; iii) the spectra were convolved for the 
instrumental and rotational broadening We used as initial 
guesses for temperatures the values derived byjPecaut, Ma-j 


majek, & Bubar (20121, that we report for the sake of clarity 
in Table [2| In Fig. we plot the synthetic profiles (even for 
Ha) superposed on the observed spectra. Moreover, we com¬ 
puted the logg and the usini by spectral synthesis of the 
region around Mgi triplet at AA5167-5183 A (see Fig. |^. 
The resulting values are listed in Table 

In a completely similar way we analysed the FEROS 
spectrum for EPIG 204506777 (see Fig. [^, which was ob¬ 
tained in better conditions with respect to the CAOS 
one (USco is high in sky in the Southern hemisphere). 
As a result we found a value of the radial velocity of 
Vrad = -13.75±0.49 km s"^ (JD = 2453906.6783). Com¬ 
parison with the value reported in Table reveals a clear 
variation due to the motion of the binary or triple system. 
Further data will be obtained in the future to build up a 
complete radial velocity curve. 

An inspection of Table reveals that there is an over¬ 
all good agreement within the errors between literature 
( Pecaut, Mamajek, fc BubarpOl^ and references therein) 
and CAOS effective temperatures. The only relevant ex¬ 
ceptions are EPIC 204399980 (HIP 79476) and 204494885 
(HIP 80130). For the first star we find an effective temper¬ 
ature about 500 K larger than that reported by [Pecaut, 


Mamajek, & Bubar (20121, who in turn refer to Vieira et 


al. (20031. However, on the basis of high-resolution spec¬ 


troscopy, ICarmoim (20101 discussed the spectral type 
of this star, concluding that it cannot be later than A5Ve and 
estimated a T®g®'^=8200 K, almost identical to our value. On 


is sound. In passing, we note that our CAOS spectrum allows 
us to evaluate the possible multiplicity of EPIC 204399980, 
by c omparing our determina tion of Vra(i=—7.4±6.4 km s“^ 
with Carmona et al. (2010l’s Vrad=0±3 km s“^. The two 


values agree within ~ 1 cr, hence we conclude that at least 
with present data, there is no evidence of multiplicity for 
this star. 

As for EPIC 204494885, the previous effective temper¬ 
ature evaluation dates back to Houk & Smith-Moore (19881 
who classified the stars visually on objective-prism plates. 
Therefore their effective temperature estimate cannot be ac¬ 
curate. Also in this case we decided to take our value. 


® Only for HIP 79476 this method has not been applied because 
of the strong emission in the Ho profile and the possible filling of 
the H/3. In this case we use the ionization equilibrium of iron to 
estimate the effective temperature 
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Figure 5. Portions of the spectral echelle orders centered around Ho and for the program stars. Superposed (red lines) the synthetic 
spectra computed with the parameters reported in Table 
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Wavelength (A) 


Figure 6. Spectral synthesis in the range of the Mgi triplet 5150 
- 5200 A 


5 EVOLUTIONARY STATUS OF THE 
INVESTIGATED STARS. 

In the context of the present paper it is important to estab¬ 
lish the evolutionary status of the target stars, particularly 
those pulsating as 5 Set or 7 Dor. Figure shows the HR 
diagram for the A- and F-type stars members of the USco 
association according to analysis of [Pecaut, Mamajek, fc| 
Bubar (20121. For the stellar parameters we have adopted 
the values listed in columns (9) and (10) of Table except 
for the four 5 Set stars for whom we used our spectroscop¬ 
ically determined values. As we discussed in Sect. there 
is some disagreement between our own results and those of 


Pecaut, Mamajek, & Bubar (20121. However, we note that 


the use of Pecaut, Mamajek, & Bubar (20121 data would 


not affect significantly our conclusions. 

The HR diagram in Figure includes some selected 
evolutionary tracks and isochrones from the models by 
Tognelli, Prada Moroni, & Degl’Innocenti (2011 1 . Note that 


Pecaut, Mamajek, & Bubar (20121 used for their analysis a 



Figure 7. Portion of the spectral echelle order in the region of 
Hcf for HIP 78977. Superposed (red line) the synthetic spectra 
computed with the parameters reported in Table 


combination of empirical isochrones and the models of|Dot-| 


ter et al. (2008 \ and found an age for USco of 11±1 Myr. As 
seen in FigTjSpwe confirm this estimate since the position of 
almost all the A-type stars is compatible within the errors 
with the 10 Myr isochrone of the models of|Tognelli, Prada] 


Moroni, & DeglTnnocenti (20111. 


So far, we have assumed that all the stars investigated 
in this paper are in the PMS phase. One can argue that 
their membership to the USco association should guarantee 
their youth, considering that not only [Pecaut, Mamajek, ^ 
|Bubar|([2012[), but also previous studies report an age of 5-10 


Myrs (see, e.g. de Zeeuw et al. 1999 Luhman & Mamajek 


2012 and references therein) for this association. However, 


as a check, we now proceed to provide further evidence on 
this assumption. 

In general, the distinctive properties of young A and F- 
type stars (also known as Herbig Ae stars) are the presence 
of Ha emission (or filling of the line) in their spectrum, a 
luminosity class HI-V, and the presence of excess infrared 


emission due to circumstellar dust (see e.g. van den Ancker, 
jde Winter, &: Tjin A Djie|1998| and reference therein). Now, 
all the stars of our sample satisfy the criterion on the lumi¬ 
nosity class. As for the presence of Ha emission, only the 
spectrum of EPIC 204399980 shows it unambiguously (cf. 
Sect. 1^ and below), whereas for the remaining stars a par¬ 
tial filling of the Ha line cannot be excluded. 

The disk population of USco has been thoroughly in- 
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vestigated by |Luhman fc Mamajek (20121. Their results for 
the stars in common with our sample are shown in column 
(11) of Table[^ About half of our stars show the presence of 
a disk, even if in most cases it is an evolved (debris, transi¬ 
tional) structure, possibly devoid of gas. The only exception 
is EPIC 204399980 with its clear signature of accretion in 
the Ha profile. Not surprisingly, its light curve displays long 
term, irregular photometric variations due to the variable 
obscuration of the photosphere produced by circumstellar 
dust, as commonly found among Herbig Ae stars (see e.g. 
Herbst et al.|1994 and references therein). An obvious inter¬ 


pretation would be that the 5 Set star EPIC 204399980 is one 
of the least evolved stars of our sample, much younger that 
the rest. However, the observed position in the HR diagram 
very close to the ZAMS runs contrary to this hypothesis, un¬ 
less the stellar luminosity has been severely underestimated 
due to a wrong extinction correction. 

Considering the other pulsating stars, according to the 
results of |Luhman fc Mamajek] ( |2012[ ), EPIC 203931628, 
204054556, and 204372172 have a debris/evolved transi¬ 
tional disk, a clear indication of their PMS nature. In the 
case of EPIC 204175508, 204494885, and 204638251 no in¬ 
frared excess was found in the literature and thus no evi¬ 
dence for the presence of a circumstellar disk. However, the 
absence of a significant infrared excess does not necessarily 
mean that these stars are evolved objects, not belonging to 
the USco association. Recent studies based both on Spitzer 
photometry up to 70 pm ( Carpenter et al.||2006 20091 and 
spectroscopy ( Dahm fc Carpenter|2009' I have shown that in 
general the disks in USco are more evolved than those ob¬ 
served in younger star forming regions. Only ~20% of the 
stars have an infrared excess at any wavelength, mostly due 
to debris rather than primordial disks. On the basis of these 
considerations, unless the membership evaluation by|Pecaut(] 


Mamajek, & Bubar (20121 is wrong, we can consider all the 


new pulsators as young intermediate-mass stars with an age 
of ~10 Myr. 


6 THE PULSATING STARS 


Recently, Zwintz et al. (20141 have observed a large sam¬ 


ple of PMS 5 Set stars and determined their effective tem¬ 
peratures and surface gravities from spectroscopic data. An 
important result of this study is that the highest observed 
p-mode frequency, Vraax, correlates quite well with the lo- 




MT^L-^ 


(see 


nil, aiaffram 

Brown et al. 

Since 

1991 

IX scales as 

Kjeldsen & 

.||2014[ for details). The coolest 


and least-evolved stars (i.e. closest to the birth-line) have 
the lowest values of i/max, while the hottest and most evolved 
stars (closest to the ZAMS) have the highest values of iZmax- 


Zwintz et al. (20141 interpret the observed trend as due to 


the behaviour of the acoustic cut-off frequency. The latter is 
directly proportional to the square root of the star’s mean 
density and therefore it should increase as a star contracts 
and its density rises. Of course, this property can only be ob¬ 
served in young stars whose radius varies with age. Indeed, 
for the more evolved, core hydrogen-burning 5 Set stars there 
is no correlation of I'max with the stellar parameters. Using 
the observed trend in a homogenous sample of pulsating 


stars in the young cluster NGC 2264, Zwintz et al. (20141 


were able to quantify the age interval during which star for¬ 
mation took place in the parent molecular cloud at the level 
of at least 5 Myr. 

As discussed in Sect. we have identihed six 
new pulsating stars: namely, EPIC 203931628, 204175508, 
204372172, 204399980, 204494885 and 204638251. Unfortu¬ 
nately, since the Nyquist frequency limit of our data is about 
24 d~^, the frequencies extracted from the LC mode data are 
subject to ambiguity preventing an accurate determination 
of t'max. On the other hand, four stars have been observed 
in SC mode and in this case we could calculate the values 
of iZmax. It turns out that in all cases iZmax ~ 42 d“^, apart 
from EPIC 204372172 with i^max ~ 33 d~^. Furthermore, all 
these stars have similar properties in terms of effective tem¬ 
peratures and luminosities. Therefore, since the relationship 
between iZmax and the stellar parameters is approximate, we 
can only draw some qualitative conclusion from the knowl¬ 
edge of the cutoff frequencies. 

First, the similar values of r'max could mean that there 
is little or no age spread in USco among the late A-type 
stars. Even if EPIC 204372172 has a smaller i/max, the posi¬ 
tion in the HR diagram is consistent with that of the other 
stars, preventing us from finding any significant difference 
in age from the comparison with the isochrones. Second, 
the fact that EPIC 204399980, the star that we consider 
to be the youngest of our sample on the basis of the pres¬ 
ence of a fully developed circumstellar disk, shows a value 
of iZmax as high as that of the other 5 Set stars may ap¬ 
pear difficult to interpret. A possibility is that the relation 


i^max-age found by Zwintz et al. 

([2014 1 is not always present 

in young clusters/associations ( 

Stabler & Palla||2014 1. Vice 


versa, and more likely, all the four pulsators of USco have 
approximately the same age, but EPIC 204399980 somehow 
managed to maintain its disk for a longer time. This is 
not an unusual situation considering that actively accreting 
stars older than 6-8 Myr and still surrounded by circum¬ 
stellar disks have been found in populous clusters, such as 


NCC 6611 ( De March! et aL]|2013 K 


In any case, this result is different from that observed 
in NCC2264 studied by ( Zwintz et al.|[^14[ ). Interestingly, 
NCC 2264 is known to have a median age of ~3 Myr, and a 


significant dispersion of the order of ~ 5 Myr (see also Dahm 


20081. In the case of USco, the age difference of its members 


seems to be very small. However, looking at their distribu¬ 
tion in the HR diagram shown in Fig. 8, we notice that there 
are several stars of USco that fall within the instability strip 
and that lie close to the 5 Myr isochrone (see the light green 
symbols). Thus, one can argue for an indication of a possible 
spread in isochronal ages in USco members of the same mag¬ 
nitude of that more solidly derived in NCC 2264. Unfortu¬ 
nately, these putative younger stars have not been observed 
by Kepler, but they are ideal targets for future dedicated 
observations to search for pulsation and therefore to probe 
the presence of a possible age spread in this association. 

It is interesting to note that low frequencies are almost 
certainly present in the PMS 5 Set stars discussed here. How¬ 
ever, one needs to be careful in making this conclusion be¬ 
cause the K2 data are not very reliable at low frequencies 
due to drift of the stars across the CCD. Nevertheless, the 
low frequencies in these stars are of sufficiently high ampli¬ 
tude that one can be confident about their identification. 
The presence of low frequencies is particularly interesting 
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Table 2. Results obtained from the spectroscopic analysis of the sample of stars presented in this work. The different columns show: 
identification, literature | |Pecaut, Mamajek, Bubar|2012[ l and spectroscopic effective temperatures, surface gravities (logg), rotational 
velocities (vsini), heliocentric Julian day of the observations and radial velocities. 


EPIC 

HIP 

(K) 

rjySpeC 

(K) 

logff 

V sin 2 
(km s“^) 

JD 

(2457100.±) 

^rad 
(km s“^) 

204506777 

78977 

6150±150 

6300 ± 200 

4.10±0.15 

80±5 

78.4706 

15.6 

± 

3.0 

204399980 

79476 

7500 ±200 

8100 ±300 

4.00±0.15 

110 ±5 

79.4551 

-7.4 

± 

6.4 

204926239 

79606 

6150±150 

6100 ± 150 

3.80±0.15 

22 ±2 

80.4255 

-24.1 

± 

1.0 

204662993 

79977 

6700±100 

6500 ± 150 

4.00±0.15 

55 ±3 

85.4380 

-6.2 

± 

2.7 

204638251 

80059 

7800 ±250 

8000 ± 250 

3.95 ±0.15 

90±5 

97.3430 

7.1 

± 

5.7 

204372172 

80088 

7500±150 

7600 ± 200 

3.80±0.15 

80±5 

96.4365 

-9.1 

± 

6.0 

204494885 

80130 

7500±150 

8000 ± 250 

3.80±0.15 

85 ±5 

78.4212 

-51.9 

± 

6.0 


in PMS stars because envelope convection is much more ex¬ 
tended in this phase, particularly at early stages when the 
star radius is still expanded. Therefore, one may expect that 
the convective blocking mechanism, responsible for the onset 
of g-modes ( Guzik et al.|20d0 l, could be more active among 
these stars. In turn, the study of these low frequency modes 
is fundamental to investigate the evolutionary status and the 
internal structure of the intermediate-mass PMS pulsators 
i Suran et al.| 2 doT l. 

The independent frequencies extracted for the 7 Dor 
star EPIC 204054556 are listed in Table and displayed 
in the bottom panel of Fig. In the same figure, we also 
show for comparison the periodgrams of the two other known 
PMS 7 Dor stars, VAS 20 and VAS 87 in the young cluster 
NGC2264 ( Zwintz et al.|20T3 K We see that the frequencies 
detected in EPIC 20405455 are systematically higher than 
the 7 Dor stars of NGC 2264, a fact that can be attributed to 
the higher effective temperature of EPIC 20405455. As can 
be seen in Fig. this star falls well within the instability 
strip for 7 Dor stars, whereas the two variables in NGC 2264 
are slightly outside (see Fig. 9 in Zwintz et al.|20i3 l. 

As a final consideration, we note that two of the in¬ 
vestigated stars, namely EPIC 204239132 and 204242194, 
do not show any pulsation at a level of ~ 0.01 ppt, albeit 
being placed inside the 5 Set instability strip, towards its 
blue boundary (see[^. The effective temperature estimate 
for these stars is not based on solid high-resolution spec¬ 
troscopy (see Pecaut, Mamajek, & Bubar 2012 and refer¬ 
ences therein). However, even a 500 K error on the effective 
temperature cannot place these stars outside the strip. The 
presence of non-variable stars inside the S Set instability 


strip is a well known long standing problem (see e.g. Balona 
|fc Dziembowski]2011a[ and references therein). A recent dis¬ 
cussion of this subject is provided by [Murphy et ah (20151 
whose investigation is based on ultra-precise Kepler pho¬ 
tometry and high-resolution spectroscopy. Their conclusion 
is that chemically normaQ non-variable stars inside the <5 Set 
instability strip exist but are rare. However, the blue limit 


of the instability strip used by these authors (i.e. Dupret 
et al.|[2005l is significantly redder than that adopted in the 


^ In A-type chemically peculiar stars, such as the Am class, a 
dumping of the pulsation is expected (see e.g. [Balona et al.|2011b[ 
and references therein). 


present work. As a result, according to [Dupret et al.| ( |2005| ) 
calculations, EPIC 204239132 and 204242194 would be close 
or outside the blue boundary of the 5 Set instability strip. 
In conclusion, given current uncertainties on both effective 
temperature measurements and instability strip determina¬ 
tions the occurrence of static stars inside the adopted pul¬ 
sation region cannot be excluded. 


7 CONCLUSIONS 

We have presented the light curves obtained with the Kepler 
spacecraft and the derived periodograms of 27 stars belong¬ 
ing to the young Upper Scorpius association. This sample 
contains only A-F type stars which are known to be subject 
to pulsation instability of the 5 Set and 7 Dor type dur¬ 
ing the contraction phase prior to the arrival on the main 
sequence. 

We have identified six new S Set variables and one 
7 Dor star (EPIC 204054556/HIP 79054) which are most 
likely PMS stars just about to arrive on the ZAMS. Four 
of the 6 Set variables, namely EPIC 204372172 (HIP 80088) 
EPIC 204399980 (HIP 79476) EPIC 204494885 (HIP 80130), 
and EPIC 204638251 (HIP 80059) were also observed in 
short cadence mode (time sampling 1 -min) allowing us to 
resolve their entire frequency spectrum. For the remain¬ 
ing two 6 Set stars, EPIC 203931628 (HIP 80196) and 
EPIC 204175508 (HIP 77960), we only got long cadence data 
(time sampling ~30-min) that were enough to unambigu¬ 
ously identify their variability. Thanks to these new data, 
USco is the second young cluster after NGC 2264 that hosts 
several PMS <5 Set /'y Dor variables discovered through space 
observations. 

One of the stars observed in short cadence, 
EPIC 204506777 (HIP 78977) appears to be a rotational 
variable in a hierarchical triple system. This is a particularly 
important system since future radial velocity measurements 
will allow us to obtain accurate determination of the mass 
of the star, thus setting strong constraints on its pulsation 
properties and on evolutionary models. 

Follow-up spectroscopic observations of the newly dis¬ 
covered variables observed in short cadence allowed us to 
obtain new estimates of Teff, logg, v sini, and Yrad- Our sec¬ 
ond epoch spectroscopy for EPIC 204506777 (the first epoch 
was obtained from archive data) reveals a significant change 
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Table 3. List of pulsation frequencies and amplitudes for the four 
candidate PMS 6 Set stars observed in SC mode. EPIC 204054556 
is a candidate 7 Dor star. Given the problematic light curve at 
low frequencies, for star EPIC 204399980, only frequencies larger 
than 10d“^ are listed. Errors on each value of frequency and 
amplitude are shown within parenthesis. 


n 

y-n 

d-1 

An 

PPt 

n 

d-1 

An 

PPt 


204638251 


204494885 

1 

26.3734(1) 

4.09(1) 

1 

22.3371(1) 

0.741(2) 

2 

22.21235) 

3.2o5) 

2 

28.68625) 

0.4455) 

3 

20.4388(1) 

2.995) 

3 

26.18065) 

0.3955) 

4 

21.25695) 

2.575) 

4 

23.11345) 

0.3325) 

5 

32.26325) 

2.575) 

5 

20.23325) 

0.2775) 

6 

27.89365) 

2.465) 

6 

32.19275) 

0.2165) 

7 

34.89665) 

2.285) 

7 

19.0010(1) 

0.19i5) 

8 

20.50555) 

2.285) 

9 

28.64595) 

0.15o5) 

9 

26.849i5) 

1.835) 

10 

21.94725) 

0.1025) 

10 

29.80555) 

1.595) 

11 

27.72895) 

0.0885) 

11 

23.16065) 

1.5o5) 

12 

18.9506(2) 

0.08i5) 

12 

19.40655) 

1.455) 

13 

29.7855 5) 

0.0795) 

13 

25.421o5) 

1.43(1) 

14 

20.38185) 

0.0795) 

14 

2 4.24545) 

1.38(1) 

15 

31.49095) 

0.070(2) 

15 

32.05985) 

1.11(1) 

16 

35.8862(3) 

0.055(2) 

16 

20.84225) 

1.06(1) 

17 

24.65525) 

0.050(2) 

17 

17.52765) 

1.045) 

18 

41.53045) 

0.0465) 

18 

27.48455) 

i.oo5) 

19 

17.80075) 

0.0455) 

19 

32.690o5) 

0.995) 

20 

28.9853(4) 

0.0445) 

20 

28.009i5) 

0.895) 

21 

35.22875) 

0.04i5) 

21 

34.31085) 

0.895) 

22 

23.53855) 

0.0375) 

22 

19.605i5) 

0.885) 

23 

25.405l5) 

0.0365) 

23 

37.96245) 

0.875) 

24 

29.9793(5) 

0.0345) 

24 

18.65225) 

0.8i5) 

25 

40.97745) 

0.0295) 

25 

19.72385) 

0.775) 

26 

31.14845) 

0.0295) 

26 

24.884i5) 

0.735) 

27 

29.45195) 

0.0295) 

27 

29.13185) 

0.685) 

28 

31.05795) 

0.027(2) 

28 

28.17135) 

0.585) 

29 

29.65885) 

0.027(2) 

29 

21.96985) 

0.585) 




30 

41.01775) 

0.565) 
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Figure 8 . HR diagram for the A-F stars members of the USco as¬ 
sociation according to [Pecaut, Mamajek, Sz Bubar| ( |2012| | except 
for the four (5 Set observed with CAOS. For these stars we used 
the results of Table Light green symbols show objects not ob¬ 
served by Kepler. Black four-starred symbols represent constant 
stars. Green triangles, red pentagons and blue circles show bi¬ 
nary/rotational, 7 Dor variables and 5 Set, respectively. For these 
variables, empty-filled and filled symbols represent data possess¬ 
ing only LG and LC-SG cadences, respectively. For comparison 
purposes we overplot the 5 Set (blue dashed lines) and the theoret¬ 
ical edges of the 7 Dor (red dotted-dashed lines) instability strips 
by |Breger k, Pamyatnykh] l |1998| | and |Warner et al.] ( |2003^ , respec¬ 
tively. We note that the mixing-length parameter adopted in the 
computation of 7 Dor instability strip was 1.87. The solid lines in 
the top panel show the ZAMS and selected evolutionary tracks 
from |Tognelli, Prada Moroni, &: Degrinnocenti] ( |2011| |. Similarly, 
in the bottom panel solid lines show isochrones for 5,10,and 15 
Myrs as taken from the same Authors. The evolutionary models 
and isochrones were calculated for a solar mixture with Z = 0.015, 
Y = 0.278 and mixing-length = 1.68. 


of center-of-mass radial velocity for this star, confirming its 
multiplicity. However, additional data are required to con¬ 
struct a complete radial velocity curve. 

We have used the newly determined stellar parameters 
of four of the J Set stars to place them in the HR diagram 
and obtain an independent estimate on their age, confirming 
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Figure 9. Top and medium panels show the schematic pe- 
riodograms of the PMS 7 Dor variables VAS 20 and VAS 87 
in NGC2244 ( |Zwintz et al.| |2013[ |. Bottom panel displays 
the periodogram of the candidate 7 Dor variable in USco, 
EPIC 204054556. The labels corresponds to the frequencies listed 
in Table The dashed lines show the instrumental frequencies. 


the value of ~10 Myr for the members of the USco associ¬ 
ation. Interestingly, their measured usini values are less or 
close to ~ 100 kms“^, thus the effect of rotation on the pul¬ 
sation frequencies should be less extreme and the frequency 
spectra more suitable for comparison with stellar pulsation 
models (e.g. Ripepi et al.||20fl| . A quantitative analysis of 
the new 5 Set and 7 Dor stars with the purpose to determine 
their masses and evolutionary stage will be the subject of a 
forthcoming paper. 
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